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In this work we systematically investigate the condensate properties, superfluid properties and 
quantum phase transitions in interacting Bose gases trapped in disordered optical potentials. We 
numerically solve the Bose-Hubbard Hamiltonian exactly for different: (a) types of disorder, (b) dis- 
order strengths, and (c) interatomic interactions. The three types of disorder studied are: quasiperi- 
odic disorder, uniform random disorder and random speckle-type disorder. We find that the Bose 
glass, as identified by Fisher et al. [Phys. Rev. B 40, 546 (1989)], contains a normal condensate 
component and we show how the three different factors listed above affect it. 



The properties of Bose systems can be strongly affected 
by interatomic interactions and disorder. In the noninter- 
acting case disorder leads to the insulating Anderson lo- 
calization [ij . If interatomic interactions are also present, 
the interplay between the interactions and disorder can 
lead to a rich and complex arrangement of different quan- 
tum phases 0. It is even possible for states to develop 
where the BEC (Bose-Einstein condensate) fraction is 
non-superfluid 

Our work has two major goals. The first is to show 
that the normal fluid BEC found by Huang and Meng [§] 
and Giorgini et al. [3] can be identified with the Bose 
glass insulating state found by Fisher et al. Q. Using 
the tight-binding Bose-Hubbard lattice model, Fisher et 
al. [5| were able to describe an insulating state in disor- 
dered lattices which exists due to the cooperative effect 
of repulsive interactions and disorder. This 'Bose glass' 
(BG) is gapless and compressible. Taking an alternative 
approach with a microscopic low-density model involv- 
ing a dilute low-temperature hard-sphere gas with ran- 
dom scatterers in the limit of weak disorder, Huang and 
Meng H and Giorgini et al. ,4] found that the superfluid 
is more depleted by the disorder than the condensate. 
That is, there exists BEC without superfluidity. 

It was shown previously using the disordered Bose- 
Hubbard model that it is possible to obtain a normal 
condensate fraction for certain parameters and that 
there is a normal condensate fraction in the Bose glass 
phase for quasiperiodic disorder Q . In this work we show 
that the parameter regime where this occurs can be iden- 
tified with the BG region described by Fisher et al. 
by showing that the BG phase and the normal conden- 
sate phase occur together as we change the interaction 
strength, disorder strength and the type of disorder. We 
also find that in the BG regime the spatial correlations go 
to zero with distance but remain significant over multiple 
lattice sites. 

Our second goal is to show how the properties of 'dirty' 
Bose gases change with: (a) the type of disorder, (b) 
the disorder strength, and (c) the interatomic interac- 
tion strength, especially with respect to the existence of 
condensate without superfluidity. We show that there are 
significant differences in the phase diagram as we change 



the type of disorder and the nature of the differences de- 
pend on the type of disorder. 

To date the majority of the work on 'dirty' bosons has 
focused on liquid ''He in porous materials such as Vycor 
glass or aerogel In recent years however, weakly- 
interacting dilute Bose gases, especially those trapped in 
optical dipole potentials, have become increasingly im- 
portant in this field [1, 1^, [13, [HI. An important ad- 
vantage of these systems is their flexibility and their 
ease of control over important experimental parameters. 
This can be seen in recent experiments in which a va- 
riety of different types of disordered optical potentials 
of controllable strength were studied e.g., quasiperiodic 
lattices speckle potentials [1, @], and lattices with 
superimposed disorder These recent experimental 

developments make an exploration of how condensate 
and superfluid properties change with the properties of 
the disorder increasingly relevant. Another advantage of 
these systems over '^He is that BEC without superfluidity 
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only occurs in very dilute systems 

The many-body quantum state for A'^ interacting 
bosons in a ID lattice with Ng lattice sites can be de- 
scribed using the ID Bose-Hubbard (BH) Hamiltonian, 



H 



■E( 

1=1 



H.c. 



U 

1=1 



1=1 



al, a 



(1) 

are the boson creation and annihilation operators 
respectively at the i-th lattice site, fii = a^Oi. J is the 
strength of tunnelling between neighboring sites and U 
the strength of the on-site interatomic interactions. The 
values of J and U can be calculated as in Ref. [13]. is 
the on-site single particle energy at site i. 

The BH Hamiltonian ^ can be solved numerically 
by diagonalization in the occupation number basis. The 
large size of this basis, D ^ {N + Ns - f )!/(A^!(iV^ - 1)!) 
restricts the iV^ and TV that can be handled using diag- 
onalization. However, diagonalization solves Eq. ([T]) in a 
numerically exact and transparent way. As will be seen, 
even in our finite systems, significant differences can be 
seen as the type of disorder, the disorder strength and 
the interaction strength change. 

Significantly, we do not make the assumption that the 
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FIG. 1: (Color online) The on-site single particle energies 
Si for: (a) The quasiperiodic lattice £i — Asin^(fcxi) for 
A = 0.3, fc = (1 + \/5)/2 and Xi = in. (b) Uniform ran- 
dom disorder for A = 0.3. (c) Speckle disorder with standard 
deviation asp; A = 2asp ~ 0.3. 30 lattice sites are shown, 
(d), (e) and (f): The probability distributions P{ei) for lat- 
tices (a), (b) and (c) respectively. The histograms show P{si) 
over 170 lattice sites and the dotted line shows the expected 
result over many realizations. 



condensate fraction, and the superfluid fraction, fs are 
the same but instead define the two separately foUowing 
the definitions used in |6|] and |5|]. fs can be calculated 
from the response of the system to a very small phase 
variation, A9{x) — Q, across the lattice. As long as 
Q is small, the increase in energy is due solely to the 
kinetic energy of the superflow. Assuming that the phase 
variation is linear, 
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in terms of the BH coefficients 0, Instead of using 
twisted boundary conditions, it is easier to use a unitary 
transformation to make a 'twisted Hamiltonian' and use 
periodic boundary conditions. 

However BEC is not defined in terms of transport prop- 
erties but as pointed out by Penrose and Onsager from 
the properties of the one-body density matrix 14|. Fol- 
lowing this, fc is taken to be the largest eigenvalue of the 
single-particle density matrix for the ground state {ipo), 
Pij = (V'o|a]ai|'0o), divided by N, the total number of 
atoms. However, because of finite size effects, the re- 
striction fc > 1/Ns exists. Due to the ID nature of the 
system this is not a true condensate. However the phys- 
ical similarity between the definition of fc here and true 
condensate and the finite size of the system means that 
it should still be useful to consider what happens in ID. 

We compare three different types of disordered lattices: 
quasiperiodic speckle [1, Q and uniform random 
It is assumed that the lattice is deep enough that disor- 
der with strength A only alters £i significantly [Hi . A = 
obviously represents a homogeneous single-periodic lat- 
tice. For uniform random disorder Si varies flatly as 



< ffi < A [see Figs.^b) and (e)]. 

The definition of A for speckle-type disorder is dif- 
ferent as there is no upper bound for e^. A speckle 
distribution can be created experimentally by scatter- 
ing laser light from a rough surface. This results in a 
random potential that obeys the probability distribution, 
P(eO = (l/(£,;))exp(-£,/(£,)) [see Figs. mc) and 

(f)] where (.) in this context denotes a mean. We model 
the speckle numerically using the method outlined in 
Refs. [1^ and Following recent experiments we de- 
fine A = 2asp, where asp is the standard deviation of 
For a speckle, (e^) — asp- It should be 
noted that by using the BH Hamiltonian ([1]), physically 
what is being modelled is not the speckle potential alone 
but the speckle plus an optical lattice. 

Quasiperiodic disorder is different from the other types 
of disorder in that it is not random but rather, com- 
pletely deterministic. It can be created by superimpos- 
ing two incommensurate lattices (that is, the ratio of the 
periodicities is irrational). Assuming that one lattice is 
much weaker than the other, the weaker lattice can be 
considered a disordered perturbation that shifts the well- 
depths of the first lattice and hence Si by at most A [see 
Fig. ma)]. Under these circumstances the quasiperiodic 
potential can be represented in the BH Hamiltonian ^ 
by Ei = Asin^ (kxi) where k is irrational and Xi = in. i 
here enumerates the lattice sites. Comparing quasiperi- 
odic disorder to random disorder is important as recent 
experiments have used quasiperiodic disorder to approx- 
imate uniform random disorder 

Fig.[2]shows how Egap, fc, fs and {fc~fs) change with: 
(1) the type of disorder, (2) the strength of the disorder 
A, and (3) the interatomic interactions U . A and U are 
both scaled by the tunnelling strength J . N = 7 and 
Ng = 7 for each sample and for the two random types 
of disorder the average of 24 samples is taken. For the 
quasiperiodic disorder we use k = {1 + ^/5)/2. 

The zero disorder line (A/J — 0) represents the ho- 
mogeneous periodic lattice. As U/J increases, the sys- 
tem undergoes a SF-MI (superfluid-Mott insulator) tran- 
sition. The MI is characterized here by: integer filling fii 
(for our parameters fii — 1), suppression of the number 
fluctuations ai, a finite energy gap, fs — and fc going 
to 1/Ns. 

In the non-interacting case {U/J ~ 0), disorder causes 
the atoms to become Anderson localized. This is char- 
acterised by the exponential spatial localization of the 
atoms and fs — 0. fc is trivially equal to one here. 

For all three types of disorder and for any A/J, as 
U/J increases the onsite interactions eventually domi- 
nate and the system becomes a MI. However, for strong 
enough disorder and moderate interactions fs is still zero 
or close to zero but the other properties of the state are 
different from that of a MI in a homogeneous lattice. 
One difference is in the fii [e.g. see Fig. |3Ib)]. However, 
the most clear difference can be seen in the energy gap 
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FIG. 2: (Color online) The energy gap Egap, the superfluid 
fraction fs, the condensate fraction /c, and (/c — fs), for the 
three different types of disordered lattices in Fig. [1] A/ J 
is the strength of the disorder and U /J is the interaction 
strength, both scaled by the tunnelling strength J . N = 7 
and Ns — 7 for each sample. For the uniform random disor- 
der and the speckle-type disorder the results shown are the 
average of 24 samples. BG='Bose glass'; MI='Mott insula- 
tor'; SF='superfluid'. 




FIG. 3: (Color online) One of the speckle samples used in 
Fig. [3 (a)-(c) A/J = 30 only, (a) The topology of the dis- 
order, Ei. The dotted line shows A/J. (b) number densities 
fii, (c) number fluctuations at, (d) the energy gap Egap, (e) 
the condensate fraction fc, (f) the superfluid fraction fa and 
(g) the difference between the condensate and superfluid frac- 
tions. 



[Figs. Ufa), (e) and (i)]. A Bose glass in a disordered 
interacting boson system is a gapless compressible insu- 
lator Q. In comparison a MI is an incompressible insu- 
lator with a gap which occurs for large enough U/J and 
commensurate filling 01 . A finite size system always has 
a discrete excitation spectrum and hence it is impossible 
to have Egap = 0. However, in this region Egap is very 
small and many times weaker than the energy gap in the 
MI which reaches values of Egap ^ 50 in Fig. [2l In the 
MI, Egap increases linearly with U/J so this is not likely 
to be its maximum value. Thus in the transition from 
the Mott insulator to this region the energy gap can be 
said to effectively disappear. Hence this region matches 
Fisher et a/.'s description of the Bose glass [2|. 

As Figs.[51Jd), (h) and (1) show, in both the superfluid 
and the MI regions there are no significant differences 
between fc and fs- However in the BG region it is clear 
that though fc is depleted it is significantly greater than 
fs. In the BG, there exists condensate without superflu- 
idity. Close to the MI or the superfluid regions, [fc — fs] 
gradually decreases. Therefore, condensate without su- 
perfluidity can be identified with the Bose glass. 



Both repulsive interactions and disorder are necessary 
for BG. However, ii U/J ^ A/J the system becomes a 
MI. As can be seen in Fig. [21 as A increases the criti- 
cal U required for the BG-MI quantum phase transition 
increases. This means that in disordered lattices using 
very strong disorder will increase the size of the param- 
eter regime where it is possible to obtain BG and hence 
condensate without superfluidity. 

The results above apply to all three types of disorder 
studied. However significant differences exist between 
the different types of disorder. First, the changes in fs 
and fc have a distinct complex lobe-like structure in the 
quasiperiodic case. This has been observed previously for 
quasiperiodic disorder and attributed to the non-random 
nature of the quasiperiodic distribution and the sensitiv- 
ity of rearrangements of the many-body wavefunction to 
the distribution of Si [^, [l3| ■ 

Another major difference exists between bounded 
(quasiperiodic and uniform random) and unbounded 
(speckle) disorder. Fig. [21 shows that for the same A/J 
using unbounded rather than bounded disorder greatly 
increases the size of the BG region. Fig. [21 shows the av- 
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FIG. 4: (Color online) The spatial correlations p{i, r) between 
sites i and i + r (with periodic boundary conditions). p{i,r = 
0) shows fii. (a)-(e) is for the speckle sample shown in Fig. [3] 
(f ) is for the quasiperiodic lattice in Fig. (2] All examples are 
for A/J = 30. 



increasing interaction strength may not occur due to the 
absence of a MI phase. An alternative is using noninteger 
N/Ns as previously studied in disordered lattices [1, [a]- 
In conclusion we have shown how superfluid and con- 
densate properties change in 'dirty' Bose gases as inter- 
atomic interactions, disorder strength and the type of dis- 
order are varied. We showed that in a disordered lattice, 
the Bose glass as identified by the lack of superfluidity 
and energy gap, comprises of a BEC without superflu- 
idity. In order to maximise the size of the parameter 
regime where it is possible to obtain BEC which is part 
of the normal fluid either the strength of the disorder can 
be increased or an unbounded disorder can be used. By 
examining the spatial correlations between lattice sites 
in the different insulating states we found that the Bose 
glass may be considered a localized homogeneous con- 
densate pinned by disorder. 



erage over many samples. Examining individual speckle 
samples shows the largest BG regions occur when there 
are "spikes" where > A. These inevitably occur due to 
the unbounded nature of the speckle distribution. An ex- 
ample of one such speckle sample with a "spike" is given 
in Fig. [3] which clearly shows the enlarged BG region. 

Studying individual samples also allows us to ex- 
amine the spatial correlations which are given by the 
off-diagonal terms of the one-body density matrix, 
Pij — {ipo\ajai\4>Q)- We let j = i + r and use peri- 
odic boundary conditions; p{i,r — 0) is thus simply 
fii. Figs, [lla)-(e) show p{i,r) for the sample disorder 
in Fig. [3l Fig. [IJf) is for the quasiperiodic disorder in 
Fig. El AU examples are for A/J = 30. 

In the non-interacting case [Fig. [D^a)] (Anderson lo- 
calization), the wavefunction is strongly localized. For 
U/J ^ A/J (Mott insulator). Hi = 1 for every lattice 
site regardless of Si [Fig. IH^e)] . No correlations exist be- 
tween lattice sites; each site is independent. 

In the BG, the density distribution is much more ho- 
mogeneous compared to the noninteracting case. Unlike 
the MI, the BG retains significant correlations between 
lattice sites. As r increases p{i,r) decreases but remains 
significant over a finite area. This also occurs in other 
disorder samples though the exact details depend greatly 
on the topology of the sample. Fig. ID^f ) shows an exam- 
ple for quasiperiodic disorder. Since fc is significantly 
greater than fs in the BG region, it is possible to con- 
sider this a localized homogeneous condensate pinned by 
the disorder. The identification with the BG hints at the 
possible structure of the quantum phase though it re- 
mains to be seen if these results apply to larger systems. 

As U/ J increases the spatial correlations between sites 
disappears. However it is possible that in a true speckle 
potential (without a lattice) the increasing fragmenta- 
tion of the condensate which can be observed here with 
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